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Extensive calculations of the electric-dipole matrix elements in alkali-metal atoms are conducted 
using the relativistic all-order method. This approach is a linearized version of the coupled-cluster 
method, which sums infinite sets of many-body perturbation theory terms. All allowed transitions 
between the lowest ns, npi/2, 'T-P3/2 states and a large number of excited states are considered in 
these calculations and their accuracy is evaluated. The resulting electric-dipole matrix elements are 
used for the high-precision calculation of frequency-dependent polarizabilities of the excited states of 
alkali-metal atoms. We find "magic" wavelengths in alkali-metal atoms for which the ns and npi/2 
and atomic levels have the same ac Stark shifts, which facilitates state-insensitive optical 

cooling and trapping. 

PACS numbers: 32.10.Dk, 32.80.Pj, 31.15.Dv, 32.70.Jz 



I. INTRODUCTION 

Recent progress in manipulation of neutral atoms in 
optical dipole traps offers advancement in a wide vari- 
ety of applications. One such application is toward the 
quantum computational scheme, which realizes qubits as 
the internal states of trapped neutral atoms [IJ . In this 
scheme, it is essential to precisely localize and control 
neutral atoms with minimum decoherence. Other ap- 
plications include the next generation of atomic clocks, 
which may attain relative uncertainty of 10~^*, enabling 
new tests of fundamental physics, more accurate mea- 
surements of fundamental constants and their time de- 
pendence, further improvement of Global Positioning 
System measurements, etc. 

In a far-detuned optical dipole trap, the potential ex- 
perienced by an atom can be either attractive or repulsive 
depending on the sign of the frequency-dependent Stark 
shift (ac Stark shift) due to the trap light. The excited 
states may experience an ac Stark shift with an opposite 
sign of the ground state Stark shift, affecting the fidelity 
of the experiments. A solution to this problem was pro- 
posed by Katori et al. 0, who suggested that the laser 
can be tuned to a magic wavelength Amagic, where lattice 
potentials of equal depth are produced for the two elec- 
tronic states of the clock transition. In their experiment, 
they demonstrated that a Amagic exists for the ^Sq—^ Pq 
clock transition of ®^Sr in an optical lattice. Four years 
later, McKeever et al. Q demonstrated state-insensitive 
trapping of Cs atoms at Amagic ~ 935 nm while still main- 
taining a strong coupling for the 6733/2 — 6si/2 transition. 
The ability to trap neutral atoms inside high-Q cavities in 
the strong coupling regime is of particular importance to 
the quantum computation and communication schemes 

i- 

In this paper, we evaluate the magic wavelengths in Na, 
K, Rb, and Cs atoms for which the ns ground state and 



either of the first two npj excited states experience the 
same optical potential for state-insensitive cooling and 
trapping. We accomplish this by matching the ac polar- 
izabilities of the atomic ns and npj states. We conduct 
extensive calculations of the relevant electric-dipole ma- 
trix elements using the relativistic all-order method and 
evaluate the uncertainties of the resulting ac polarizabil- 
ities. We also study the ac Stark shifts of these atoms 
to determine the dependence of Amagic on their hyperfine 
structure. 

The paper is organized as follows. In section[TTl we give 
a short description of the method used for the calculation 
of the ac polarizabilities and list our results for scalar and 
tensor ac polarizabilities. In section ITTTl we discuss the 
effect of ac Stark shifts on the hyperfine structure of the 
alkali- metal atoms. In section HVl we discuss the magic 
wavelength results for each of the atoms considered in 
this work. 



II. DYNAMIC POLARIZABILITIES 

We begin with an outline of calculations of the ac 
Stark shift for linearly polarized light, following Refs. [1, 
Q. Angels and Sandars 3 discussed the methodology 
for the calculation of the Stark shift and parameteriza- 
tion of the Stark shift in terms of the scalar and tensor 
polarizabilities. Stark shifts are obtained as the energy 
eigenvalues of the Schrodinger equation with interaction 
operator V/ given by 



Vi 



-e • d. 



(1) 



where e is the applied external electric field and d is the 
electric-dipole operator. The first-order shift associated 
with Vi vanishes in alkali-metal atoms. Therefore, the 
Stark shift Ai? of level v is calculated from the second- 
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order expression 



k 



E„ — El. 



(2) 



where the sum over k includes all intermediate states 
allowed by electric-dipole transition selection rules, and 
Ek is the energy of the state k. 

Using the Wigner-Eckart theorem, one finds that AE 
can be written as the sum M 



AE 



- + 1) 
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where ao{uj) and 02(0;) are the scalar and tensor ac po- 
larizabilities, respectively, of an atomic state v. The 
laser frequency uj is assumed to be several lincwidths 
off-resonance. Here, the polarization vector of the light 
defines the z direction. 

The scalar ac polarizability ao(w) of an atom can be 
further separated into an ionic core contribution acorc('-^) 
and a valence contribution aQ(w). The core contribu- 
tion has a weak dependence on the frequency for the 
values of oj relevant to this work. Therefore, we use the 
static ionic core polarizability values calculated using the 
random-phase approximation (RPA) in Ref. Q. The va- 
lence contribution ao(^) static polarizability of a 
monovalent atom in a state v is given by Q 



3(2j, + 1) 



E 



{k\\d\\vy{Ek^E, 
[Ek-E,Y-uj^ 



(4) 



where (fcHdHw) is the reduced electric-dipole (El) ma- 
trix element. The experimental energies Ei of the most 
important states i which contribute to this sum have 
been compiled for the alkali atoms in Refs. [1, [l^, 
Unless stated otherwise, we use atomic units (a.u.) for 
all matrix elements and polarizabilities throughout this 
paper: the numerical values of the elementary charge, 
e, the reduced Planck constant, fi — /i/27r, and the 
electron mass, me, are set equal to 1. The atomic 
unit for polarizability can be converted to SI units via 
a/h [Hz/(V/m)2]=2.48832xl0-*a [a.u.], where the con- 
version coefficient is 47reoao/^ and the Planck constant h 
is factored out. 

The tensor ac polarizability 0:2 (w) is given by p^ : 



a2(w) 



C 
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' (Ek-E^r-Lu^ ' 

5j.(2j. -1) 



(5) 



6(j. + l)(2j„ + l)(2j,+3) 
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The ground state ac polarizabilities of alkali-metal 
atoms have been calculated to high precision [l^. How- 
ever, no accurate systematic study of the ac polarizabili- 
ties of the excited states of alkali-metal atoms is currently 



TABLE I: Contributions to the 5p3/2 scalar (ao) and ten- 
sor (02) polarizabilities at A=790 nm in Rb and their uncer- 
tainties in units of Oq. The absolute values of corresponding 
reduced electric-dipole matrix elements d (in eao) and the cor- 
responding transition wavelength in vacuum Arcs (in nm) are 
also given. 
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Ct2 


5p3/2 


— i3Si/2 


( oU 


o.y ( ( 


-4i0o(^0 ) 


4i0o(^0 


5p3/2 


— GSi/2 


1367 


6.047 


-92(1) 


92(1 


5p3/2 


1/ z 


741 


1.350 


41.1(1) 


-41.1(1 


5p3/2 


— 8Si /2 


616 


0.708 


2.88(2) 


-2.88(2 


5P3/2 


— 9Sl /2 


566 


0.466 


0.922(6) 


-0.922(6 
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— 10s, /2 
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-0.430(3 


Op-i/2 


A J 

- 4d3/2 


1 con 
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— 5^3/2 
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— 6C/3/2 
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0.506 
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1.3(3 
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- 7^3/2 
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0.370 
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0.48(7 


5p3/2 


— 8(^3/2 
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0.283 


0.30(4) 


0.24(3 
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0.14(2 
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- 4d5/2 
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10.899 


-242(4) 


48.4(8 


5p3/2 


- 5^5/2 


776 


1.983 


317(28) 


-63(6 


5p3/2 


- 6^5/2 


630 


1.512 


14(3) 


-2.9(6 


5p3/2 


- 7^5/2 


573 


1.104 


5.4(8) 


-1.1(2 


5p3/2 


— 8d5/2 


543 


0.845 


2.7(3) 


-0.54(6 


5p3/2 


- 9^5/2 


526 


0.672 


1.6(2) 


-0.31(3 


Otail 








19(14) 


-5(5 


dcorc 








9.1(5) 




Total 








-4060(32) 


4184(9 



available. The polarizability calculations for the excited 
np states are relatively complicated because in addition 
to p — s transitions they also involve p — d transition 
matrix elements. The matrix elements involving the nd 
states are generally more difficult to evaluate accurately, 
especially for the heavier alkalies. 

In this work, we calculate np — n'd transition matrix 
elements using the relativistic all-order method [3, 
and use these values to accurately determine the ripi/2 
and np^/2 state ac polarizabilities. In the relativistic all- 
order method, all single and double (SD) excitations of 
the Dirac-Fock (DF) wave function are included to all or- 
ders of perturbation theory. For some matrix elements, 
we found it necessary to also include single, double and 
partial triple (SDpT) excitations into the wave functions 
(SDpT method). We conduct additional semi-empirical 
scaling of our all-order SD and SDpT values where we 
expect scaled values to be more accurate or for more ac- 
curate evaluation of the uncertainties. The scaling pro- 
cedure has been described in Refs. [sl. ITsl [l6t . 

We start the calculation of the np state valence polar- 
izabilities using Eqs. ^ and ([5]). For the wavelength 
range considered in this work, the first few terms in the 
sums over k give the dominant contributions. Therefore, 
we can separate the np state valence polarizability into 
a main part, amain, that includes these dominant terms, 
and a remainder, ataii- We use a complete set of DF 
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wave functions on a nonlinear grid generated using B- 
splines [13] in all our calculations. We use 70 splines of 
order 11 for each value of the angular momentum. A 
cavity radius of 220 a.u. is chosen to accommodate all 
valence orbitals of amain- In our K and Rb calculations, 
we include all ns states up to 10s and all nd states up to 
9d; lis, 12s, and 10c? are also added for Cs. Such a large 
number of states is needed to reduce uncertainties in the 
remainder ataii- We use the experimental values com- 
piled in Ref. [la| along with their uncertainties for the 
first np — ns matrix elements, for example the 5pj — 5s 
matrix elements in Rb. We use the SD scaled values for 
some of the np — n'd and np — n's matrix elements in the 
cases where it was essential to reduce the uncertainty 
of our calculations and where the scaling is expected to 
produce more accurate results based on the type of the 
dominant correlation corrections. This issue is discussed 
in detail in Refs. [l^, [l^] and references therein. 

In Table [H we give the contributions to the scalar and 
tensor polarizabilities of the Rb 5p3/2 state at 790 nm 
to illustrate the details of the calculation. The absolute 
values of the corresponding reduced electric-dipole ma- 
trix elements, d, used in the calculations are also given. 
The contributions from the main term are listed sepa- 
rately. We also list the resonant wavelengths Arcs corre- 
sponding to each transition to illustrate which transitions 
are close to 790 nm. As noted above, we use the exper- 
imental values for the 5p3/2 — 5s matrix element from 
the Ref. [31 . We use the recommended values for the 
5p3/2 ~ 4(ij transitions derived from the Stark shift mea- 
surements [111 in Ref. [13] . We find that the contribution 
of the 5p3/2 — 5s transition is dominant since the wave- 
length of this transition (Arcs = 780 nm) is the closest 
to the laser wavelength. The next dominant contribu- 
tion for the scalar polarizability is from the 5p3/2 — 5c?5/2 
transition (Aros = 776 nm). While the contribution from 
this transition is less than one tenth of the dominant 
contributions, it gives the dominant contribution to the 
final uncertainty owing to a very large correlation cor- 
rection to the 5p3/2 ~ 5(i5/2 reduced electric-dipole ma- 
trix element. In fact, the lowest-order DF value for this 
transition is only 0.493 a.u. while our final (SD scaled) 
value is 1.983 a.u. We take the uncertainty in this tran- 
sition to be the maximum difference of our final values 
and ab initio SDpT and scaled SDpT values. While the 
5p3/2 ~ transition has almost the same transition 

wavelength owing to the very small fine-structure split- 
ting of the 5d state, the corresponding contribution is 
nine times smaller owing to the fact that the 5^3/2 — 5^3/2 
reduced electric-dipole matrix element is smaller than the 
5p3/2 — 5(^5/2 matrix element by a factor of three. As ex- 
pected, the contributions from the core and tail terms 
are very small in comparison with the total polarizabil- 
ity values at this wavelength. 

In Table ini we compare our results for the first excited 
npi/2 and np^/2 state static polarizabilities for Na, K, 
Rb, and Cs with the previous experimental and theoreti- 
cal studies. The measurements of the ground state static 



TABLE II: Comparison of static polarizabilities of npi/2 and 
np3/2 states with other experiments and theory. "Ref. [23l |. 
''Ref. [Ill, '^Ref. ''derived from Ref. [H] Dl line Stark 
shift measurements and recommended values for ground state 
polarizability from Ref. [13, ''Ref. [H, ■''Ref. [Q, "derived 
from Ref. [2l[ Dl line Stark shift measurement and ground 
state polarizability measurement from Ref. [30|, ''derived from 
Ref. 31] D2 line Stark shift measurement and ground state 
polarizability from Ref. [s^]- Units: ag. 



Na 


ao(3pi/2) 


ao(3p3/2) 


Q!2(3p3/2j 


Present 


359.9 


361.6 


-88.4 


Other 


359.7" 


361.4" 


-88.0" 


Exp. 


359.2(6)'' 


360.4(7)'' 


-88.3(4)'= 


K 


ao(4pi/2) 


ao(4p3/2) 


Q!2(4p3/2) 


Present 


602 


613 


-109 


Other 


605" 


616" 


-111" 


Exp. 


606.7(6)'' 


614(10)'= 


-107(2)'= 


Rb 


ao(5pi/2) 


ao(5p3/2) 


Q!2(5p3/2) 


Present 


805 


867 


-167 


Other 


807" 


870" 


-171" 


Exp. 


810.6(6)'' 


857(10)'= 


-163(3)'= 


Cs 


ao(6pi/2) 


ao(6p3/2) 


Q2(6p3/2) 


Present 


1338 


1650 


-261 


Other 


1290-'' 


1600^ 


-233^ 


Exp. 


1328.4(6)3 


1641(2)'' 


-262(2)'' 



polarizability of Na by Ekstrom et al. [241] were combined 
with the experimental Stark shifts from Refs. [1^, [3^ to 
predict precise values for the 3pi/2 and 3p^/2 scalar po- 
larizabilities [13] ■ The tensor polarizability of the 3p3/2 
state of Na has been measured by Windholz et al. [2^ . 
The Stark shift measurements for K and Rb have been 
carried out by Miller et al. [1^ for Dl lines and by Krenn 
et al. [IH] for D2 lines. We have combined these Stark 
shift measurements with the recommended ground state 
polarizability values from Ref. (27j to obtain the npj po- 
larizability values that we quote as experimental results. 
The np^/2 tensor polarizabilities in K and Rb were mea- 
sured in Ref. [1^. Accurate Dl and D2 Stark shift mea- 
surements for Cs have been reported in Refs. [U, |3l[. 
The most accurate experimental measurement of the 6s 
ground state polarizability from Ref. [13] has been used 
to derive the values of the &P1/2 and 6^3/2 state polar- 
izabilities in Cs quoted in Table |TT1 Our results are in 
excellent agreement with the experimental values. 

We note that we use our theoretical values for the 
Ap — 3d transitions in K and 5p — Ad transitions in Rb 
to establish the accuracy of our approach. We use more 
accurate recommended values for these transitions de- 
rived from the experimental Stark shifts [ll| in Ref. [12] 
in all other calculations in this work, as described in the 
discussion of Table [J 
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III. AC STARK EFFECT FOR HYPERFINE 
LEVELS 



hyperfine level of ground ns state. 



In the above discussion, we neglected the hyperfine 
structure of the atomic levels. However, it is essential 
to include the hyperfine structure which is affected by 
the presence of external electric field for practical appli- 
cations discussed in this work. In this section, we calcu- 
late the eigenvalues of the Hamiltonian H representing 
the combined effect of Stark and hyperfine interactions. 
Then, we subtract the hyperfine splitting from the above 
eigenvalues to get the ac Stark shift of a hyperfine level. 
This value is used to calculate the ac Stark shift of the 
transition from a hyperfine level of excited np state to a 



J 



A. Matrix elements of the Stark operator 

First, we evaluate the matrix elements of Stark opera- 
tor in the hyperfine basis. The energy difference between 
two hyperfine levels is relatively small for cases consid- 
ered in this work, and the hyperfine levels are expected 
to mix even if small electric fields are applied. Therefore, 
the Stark operator now has non-zero off-diagonal matrix 
elements. The general equation for the matrix elements 
is given by 



/^Ef f" = 



k 



{IjFM \Vi\IkjkFkMk) {IkJkFkMk \Vi\ I" f F" M") 
E„ — Eh 



(6) 



where \IjFM) represents atomic states in hyperfine ba- 
sis, / is the nuclear spin, and F = I + j. 

The interaction operator V/ given by Eq. ^ does not 
affect the nuclear spin /. In addition, the shifts due to Vi 
are not high enough to cause mixing between two levels 
with different angular momentum j . As a result, AEf,f" 
is diagonal in / and j. These approximations enable us 
to label the states in hyperfine basis as \FM), and Eq. 
^ can be simplified as 



AEf^f" = J2 

k 



{EM \Vi\ FkMk) {FkMk \Vi\ F" M") 



(7) 



Ey — Ek 

One can write the above matrix element as 

AEf.f" = {FM \Vh\ F"M") , (8) 



where the Stark shift operator V// is defined in terms of 
A operator as 



Vn = ViWi, 

\FkMk){FkMk\ 



A 



Ev — Ek 



(9) 
(10) 



If the applied electric field is in the z direction, then 
the energy shifts are diagonal in AI. Thus, the matrix 
elements can be written as 



AEf^f" = {FM\Vii\F"M) 



(11) 



We use the Wigner-Eckart theorem to carry out the 
angular reduction, i.e. sum over the magnetic quantum 
numbers. Then, the matrix elements can be written in 



terms of scalar and tensor polarizabilities as 



{FM\Vii\F"M) 



ao{uj)e^SF,F" 



-a2{iu)e^ {FM\Q\F"M). {12) 



The first term in the equation above, containing the 
scalar polarizability, results in the equal shifts of all of 
the hyperfine levels and is non-zero only for the diagonal 
matrix elements {F = F"). The tensor part mixes states 
of different E through Q operator. The non-zero matrix 
elements of Q are 



(j + l)(2j + l)(2j + 3) 



{EM\Q\F"M)^ 

V j(2j-l) 

x{-iy+^+^-^"-^"^{2E+l){2E" 

E 2 E" 

M -M 



1/2 



1)) 





2 






/ 


7} 



(13) 



For each magnetic sublevel, there is a matrix with rows 
and columns labeled by F and E" . Therefore, magnetic 
sublevels with different values of \M\ are shifted by a 
different amount. A detailed discussion of this matrix is 
given by Schmieder 0]. 



B. Energy eigenvalues 

Since the Stark interactions considered in this work 
are comparable to the hyperfine interactions, we find the 
combined shift of a hyperfine level by diagonalizing the 
Hamiltonian given by 



(14) 
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where Vms is the hyperfine interaction operator. In the 
hyperfine basis, l^ifs is diagonal with the foUowing matrix 
elements jssj 

{FM\Vms\FM) = ^Az + 

3z{z+l)-4I{I + l)j{j + l) 

4J(2/-l)2j(2j-l) ' ^ 

where z = f(-F+l) -/(/+!) -.7(7 + 1), and A and B are 
hyperfine-structure constants [3J]. The matrix elements 
of H which describe the combined effect of the Stark 
interaction V77 and hyperfine interaction Vhfs are given 
by 

Vf,f"-m = {FM\Vii\F"M) 

+ {FM\Vms{F ^ F")\F"M) . (16) 

Using Eq. IT^ and Eq. (fT5)) . the above matrix elements 
can be reduced to a more useful form 

Vf,F";M — —-^Cioe^SF,F" 

- ^a2e^ {FM\Q\F"M) 

+ {FM\Vms{F^F")\F"M). (17) 

The combined shift of a hyperfine level is evaluated 
by diagonalizing the matrix formed with Vf,f";M- The 
resulting diagonal matrix element {/S.Ef,f) corresponds 
to the shift in a hyperfine level -F, resulting from two 
effects: the hfs interaction T4ifs and the Stark effect Vn- 
Consequently, we should subtract the hyperfine splitting 
from the these shifts to get the ac Stark shift of a level 
given by 

Ani.FM = ^Ef^f - {FM \Vms{F = F")\ F"M) . (18) 

The ac Stark shift of the transition from an excited 
state to the ground state AE{n'l'j,F'M' nljFM) is 
determined as the difference between the ac Stark shifts 
of the two states. We calculate the magic wavelength 
where the ac Stark shift of the np — ns transition is equal 
to zero. The results of the calculation are presented in 
the next section. 



IV. MAGIC WAVELENGTHS FOR THE np - ns 
TRANSITIONS 



We define the magic wavelength Amagic as the wave- 
length where the ac polarizabilities of the two states are 
the same, leading to zero ac Stark shift for a correspond- 
ing transition. For np — ns transitions considered in this 
work, it is found at the crossing of the ac polarizabil- 
ity curves for the ns and np states. In the case of the 
np^/2 ~ n-s transitions, the magic wavelengths need to 
be determined separately for the cases with mj = ±1/2 



TABLE III; Magic wavelengths Amagic above 500 nm for the 
3pi/2 — 3s and 3p3/2 — 3s transition in Na and the correspond- 
ing values of polarizabilities at the magic wavelengths. The 
resonant wavelengths Arcs for transitions contributing to the 
3pj ac polarizabilities and the corresponding absolute values 
of the electric-dipole matrix elements are also listed. The 
wavelengths (in vacuum) are given in nm and electric-dipole 
matrix elements and polarizabilities are given in atomic units. 



Transition: 3pi/2 ^ 3s 



Resonances 


d 


'^rcs 




(-^magic ) 


3pi/2 


- 4s 


3.576(1) 


1138.46 


1028.7(2) 


241(1) 


3pi/2 


— 3^3/2 


6.791 


818.55 










615.88(1) 


1909(2) 


3pi/2 


- 5s 


0.757 


615.59 






3Pl/2 


-3s 


3.525(2) 


589.76 












589.457 


52760(100) 


3pi/2 


- 4^3/2 


1.917 


568.42 










566.57(1) 


-1956(3) 


3pi/2 


-6s 


0.391 


515.03 












514.72(1) 


-514(1) 



Transition: 3p3/2'^j — 3s |m 



3p3/2 


- 4s 


5.067(1) 


1140.69 


1/2 


984.8(1) 


252(1) 


3p3/2 


- 3^5/2 


9.122(1) 


819.71 








3p3/2 


— 3(^3/2 


3.041 


819.70 


1/2 


616.712(1) 


1854(2) 


3p3/2 


- 5s 


1.071 


616.25 














1/2 
3/2 


589.636 
589.557(1) 


-66230(80) 
-42(2) 


3p3/2 


-3s 


4.984(3) 


589.16 








3p3/2 


- 4^3/2 


0.857 


568.98 








3p3/2 


- 4d5/2 


2.571 


568.98 












1/2 
3/2 


567.43(1) 
566.79(1) 


-2038(3) 
-1976(3) 


3p3/2 


-6s 


0.553 


515.48 














1/2 


515.01(1) 


-517(1) 



and nij — ±3/2 owing to the presence of the tensor con- 
tribution to the total polarizability of np^^2 state. Ac- 
cording to Eq. ([3]), the total polarizability for the np^/2 
states is determined as a = ao — a2 for nij ~ ±1/2 and 
a = ao + a2 for nij — ±3/2. The uncertainties in the 
values of magic wavelengths are found as the maximum 
differences between the central value and the crossings of 
the ttns ± Sans and Unp ± Sanp curves, where the da are 
the uncertainties in the corresponding ns and np polar- 
izability values. 

We also study Amagic for transitions between particu- 
lar np3/2-F"M' and nsFM hyperfine sublevels. The ac 
Stark shifts of the hyperfine sublevels of an atomic state 
are calculated using the method described in the previ- 
ous section. In alkali-metal atoms, all magnetic sublevels 
have to be considered separately; therefore, a Amagic is 
different for the np F'M' — ns FM transitions. We in- 
clude several examples of such calculations. 

We calculated the Amagic values for npi/2 — ns and 
np3/2 — ns transitions for all alkali atoms from Na to Cs. 
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As a general rule, we do not list the magic wavelengths 
which are extremely close to the resonances. Below, we 
discuss the calculation of the magic wavelengths sepa- 
rately for each atom. The figures are presented only for 
states as they are of more experimental relevance. 
All wavelengths are given in vacuum. 



A. Na 

We list the magic wavelengths Amagic above 500 nm for 
the 3pi/2~3s and 3p3/2~3s transitions in Na and the cor- 
responding values of polarizabilities at the magic wave- 
lengths in Table IIIII For convenience of presentation, 
we also list the resonant wavelengths Arcs for transitions 
contributing to the 3pi/2 and 3p3/2 ac polarizabilities and 
the corresponding values of the electric-dipole matrix el- 
ements along with their uncertainties. Only two transi- 
tions contributing to the ground state polarizabilities are 
above 500 nm, 3pi/2 — 3s and 3p3/2 —3s. Therefore, there 
is no need to separately list resonant contributions to the 
ground state polarizability. To indicate the placement of 
the magic wavelength, we order the lists of the resonant 
and magic wavelengths to indicate their respective place- 
ment. The polarizabilities and their uncertainties are 
calculated as described in Section |TT1 The transitions up 
to np — 3s and 3p — nd, n = 6 are included into the main 
term and the remainder is evaluated in the DF approx- 
imation. The values of the 3p — 3s matrix elements are 
taken from [l8l | , the remaining matrix elements are either 
SD or SD scaled values. The uncertainties in the values 
of the Na matrix elements were estimated to be generally 
very small. The resonant wavelength values are obtained 
from energy levels from National Institute of Standards 
and Technology (NIST) database ,9,]. We assume no un- 
certainties in the energy values for all elements. 

Since the 3s polarizability has only two resonant tran- 
sitions at wavelengths greater than 500 nm, it is generally 
small except in close vicinity to those resonances. Since 
the polarizability of the 3pi/2 state has several contri- 
butions from the resonant transitions in this range, it is 
generally expected that it crosses the 3s polarizability 
in between of the each pair of resonances listed in Ta- 
ble lllll unless the wavelength is close to 3p— 3s resonances. 
The same is expected in the case of the 3p3/2 polariz- 
ability for the \mj\ = 1/2 cases as described by Eq. ^ 
{a = ao — a2)- However, when a — ag +a2 {rrij = ±3/2 
in Eq. [3|) all 3p3/2 — ns transitions do not contribute to 
the total polarizability owing to the exact cancellation 
of the scalar and tensor contributions for v = 'ip^/2 and 
fc = ns in Eqs. (|4]) and ([5]). In this case, the angular fac- 
tor in Eq. ([5]) is exactly — 2/3(2j„ -I- 1) leading to exact 
cancellation of such terms, and the total polarizability 
comes from the remaining 3p — nd contributions which 
do not cancel out. As a result, there are no resonances for 
the TTij = ±3/2 cases at the wavelengths corresponding 
to 3p3/2 — ns transitions leading to substantial reduc- 
tion in the number of magic wavelengths. We note that 



there is a magic wavelength at the 589.557 nm owing to 
the resonances in the ground state polarizability. The 
corresponding polarizability value is very small making 
this case of limited practical use. While there has to be a 
magic wavelength between 3^3/2-4^3/2 and 3p3/2 — 4c?5/2 
resonances at 568.98 nm, we are not listing it owing to 
a very small value of the 4d fine-structure splitting. We 
illustrate the magic wavelengths for the 3p3/2 — 3s tran- 
sition in Fig. [1] where we plot the values of the ac polar- 
izabilities for the ground and 3p3/2 states. 

It is interesting to consider in more detail the region 
close to the ipj — 3s resonances since in this case one 
magic wavelength is missing for both 3pi/2 — 3s and 
3p3/2 — 3s ao — «2 cases, one on the side of each 3p — 3s 
resonance as evident from Table IIIII We plot the ac 
polarizabilities for the 3s, 3pi/2, and 3p3/2TOj states in 
this region in Fig. [2l The placements of the 3pi/2 — 3s 
and 3p3/2 — 3s resonances are shown by vertical lines. 
In the case of the 3pi/2 state, the 3pi/2 — 3s resonance 
contributes to both ground state and 3pi/2 polarizabili- 
ties. As a result, both of these polarizabilities are large 
but have opposite sign right of the ipi/2 — 3s resonance 
at 589.76 nm leading to missing magic wavelength for 
3pi/2 — 3s transition between the 3pi/2 — 5s and 3pi/2 — 3s 
resonances. In the 3p3/2 — 3s ao — ci2 case, there is a 
missing magic wavelength to the left of the 3p3/2 — 3s 
589.12 nm resonance for the same reason. The values of 
the ao + oi2 for the 3p3/2 state are very small and neg- 
ative in that entire region owing to the cancellations of 
the 3p — 3s contributions in the scalar and tensor 3p3/2 
polarizabilities described above. 

In summary, there is only one case for Na in the con- 
sidered range of the wavelengths where the magic wave- 
length exists for all sublevels (567nm) at close values of 
the polarizabilities (-2000 a.u.) The ac polarizabilities for 
the 3s and 3p3/2 states near this magic wavelength are 
plotted in Fig. [31 The plot of the ac Stark shift for the 
transition between the hyperfine sublevels near 567 nm is 
shown in Fig. 21 The Amagic is found at the point where 
the ac Stark shift of the transition from 3^3/2^^' — 3M' 
sub levels to isFM sub levels crosses zero. This cross- 
ing of ac Stark shift curve occurs close to 567 nm which 
is close to the wavelength predicted by the crossing of 
polarizabihties illustrated by Fig. [31 as expected. 



B. K 

The magic wavelengths Amagic above 600 nm for the 
Api/2 — 4s and 4p3/2 — 4s transitions in K are listed in 
Table [TVl Table [Tvl is structured in exactly the same way 
as Table IIIII The electric-dipole matrix elements for the 
4p — 4s transitions are taken from [l^ , and the electric- 
dipole matrix elements for the 4p — 3d are the recom- 
mended values from Ref. |22|| derived from the accurate 
Stark shift measurements [26[ . The resonant wavelengths 
are obtained from the energy levels compiled in the NIST 
database (9|]. The transitions up to 4p — 10s and 4p — 9c? 
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FIG. 1: Frequency-dependent polarizabilities of Na atom in the ground and 3p3/2 states. The arrows show the magic wave- 
lengths. 
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FIG. 2: Magic wavelengths for the 3pi/2 — 3s and 3p3/2 — 3s 
transition of Na. 
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FIG. 3: Magic wavelength for 3p3/2 — 3s transition of Na. 



are included into the main term of the polarizabiUty, and 
the remainder is evaluated in the DF approximation. In 
the case of some higher states, such as 9s, we did not 
evaluate the uncertainties of the matrix elements where 
we expect them to be small (below 0.5%). As a result, 
the uncertainties in the values of the magic wavelengths 
near these transitions do not include these contributions 
and may be slightly larger than estimated. In the test 
case of Rb, the uncertainties are evaluated for all tran- 
sitions with resonant wavelengths above 600 nm and no 
significant differences in the uncertainties of the relevant 
magic wavelengths with other elements are observed. 

The main difference between the Na and K calculation 
is extremely large correlation correction to the values of 
the Ap — 4d transitions. The correlation correction nearly 
exactly cancels the lowest-order DF value leading to a 
value that is essentially zero within the accuracy of this 
calculation. As a result, we do not quote the values for 



the magic wavelength between Ap — Ad and 4p — 6s res- 
onances. We note that these two resonances are very 
closely spaced (1.5 nm), thus probably making the use of 
such a magic wavelength impractical. Our present cal- 
culation places the magic wavelength for the Api/2 — 4s 
transition in the direct vicinity(within 0.01 nm) of the 
693.82 nm resonance. We note that the measurement 
of the ac Stark shift (or the ratio of the 4s to 4p Stark 
shifts) near the Ap — Ad resonance may provide an excel- 
lent benchmark test of atomic theory. This problem of 
the cancellation of the lowest and higher-order terms for 
the np — nd transitions is unique to K. In the case of Rb, 
the correlation for the similar 5p — 5d transition is very 
large but adds coherently to the DF values. As a result, 
we were able to evaluate the corresponding Rb 5p — 5d 
matrix elements with 4.5% accuracy. The accuracy is 
further improved for the 6p — 6d transitions in Cs. 

We also located the magic wavelengths for the 4p3/2 — 
4s transition between 4p3/2 — 3c?3/2 and 4p3/2 — 3(i5/2 
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FIG. 4: ac Stark shifts for the transition from 3p3/2J"' = 3M' 
sub levels to 3sFM sub levels in Na as a function of wave- 
length. The electric field intensity is taken to be 1 MW/cm'^. 



TABLE IV: Magic wavelengths Amagic above 600 nm for the 
4pi/2 — 4s and 4p3/2 — 4s transition in K and the correspond- 
ing values of polarizabilities at the magic wavelengths. The 
resonant wavelengths Arcs for transitions contributing to the 
Apj ac polarizabilities and the corresponding absolute values 
of the electric-dipole matrix elements are also listed. The 
wavelengths (in vacuum) are given in nm and electric-dipole 
matrix elements and polarizabilities are given in atomic units. 

Transition: 4pi/2 — 4s 



Resonances 


d 


'^rcs 




-^niagic 


Q^(Amagic ) 


4pi/2 


- 5s 


3.885(19) 


1243.57 
















1/2 


1227.7(2) 


472(1) 


4pi/2 


— 3^3/2 


7.984(35) 


1169.34 








4Pl/2 


- 4s 


4.102(5) 


770.11 














1/2 


768.413(4) 


20990(80) 


4pi/2 


- 4d3/2 


0.097(57) 


693.82 








4pi/2 


- 6s 


0.903 


691.30 














1/2 


690.15(1) 


-1186(2) 



Transition: 4p3/2n^j — 4s 



resonances, but found that nij = ±1/2 curve crosses the 
4s polarizability very close (within 0.002 nm) to the reso- 
nance. Therefore, we do not Hst this crossing in Table HVl 
We note that rrij = ±3/2 curve crosses the 4s polariz- 
ability curve further away from resonance at 1177.35 nm. 
The polarizability values for both of these crossings is 
500 a.u. 



TABLE V: Magic wavelengths Amagic above 600 nm for the 
5pi/2 — 5s transition in Rb and 6pi/2 — 6s transitions in Cs 
and the corresponding values of polarizabilities at the magic 
wavelengths. The resonant wavelengths Arcs for transitions 
contributing to the npj ac polarizabilities and the correspond- 
ing absolute values of the electric-dipole matrix elements are 
also listed. The wavelengths (in vacuum) are given in nm and 
electric-dipole matrix elements and polarizabilities are given 
in atomic units. 

Rb Transition: 5pi/2 — 5s 



Resonances 


d 


Arcs 


magic 


'^('^magic) 


5pi/2 


- 4d3/2 


8.051(67) 


1475.65 














1350.9(5) 


476(1) 


5pi/2 


— 6s 


4.146(27) 


1323.88 






5pi/2 


- 5s 


4.231(3) 


794.98 












787.6(1) 


5417(25) 


5pi/2 


- 5d3/2 


1.35(7) 


762.10 










761.5(1) 


-5230(30) 


5pi/2 


- 7s 


0.953(2) 


728.20 














727.35(1) 


-1877(3) 


5pi/2 


- 6d3/2 


1.07(11) 


620.80 














617.7(7) 


-490(3) 


5pi/2 


- 8s 


0.502(2) 


607.24 














606.2(1) 


-444(1) 


Cs Transition: 6p 


1/2 - 6s 








6pi/2 


- 5d3/2 


7.016(24) 


3011.15 














1520(3) 


583(2) 


6pi/2 


- 7s 


4.236(21) 


1359.20 






6pi/2 


- 6s 


4.489(7) 


894.59 






6pi/2 


- 6d3/2 


4.25(11) 


876.38 






6Pl/2 


- 8s 


1.026 


761.10 












759.40(3) 


-1282(3) 


6pi/2 


- 7d3/2 


2.05(2) 


672.51 














660.1(6) 


-513(3) 


6pi/2 


- 9s 


0.548 


635.63 












634.3(2) 


-424(2) 


6pi/2 


~ 8d3/2 


1.30(2) 


601.22 







4p3/2 


- 5s 


5.535(26) 1252.56 


1/2 


1227.7(2) 


472(1) 




4p3/2 


- 3d5/2 


10.741(47) 1177.61 










4P3/2 


— 3d3/2 


3.580(16) 1177.29 


1/2 


769.432(2) 


-27190(60) 










3/2 


768.980(3) 


-356(8) 
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4p3/2 
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5.800(8) 766.70 
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a' 
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4p3/2 
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0.033(47) 696.61 
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FIG. 5: Magic wavelengths for the 5p3/2 — 5s transition of 
Rb. 
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TABLE VI: Magic wavelengths Amagic above 600 nm for the 5p3/2 — 5s transition in Rb and 6^3/2 — 6s transitions in Cs and the 
corresponding values of polarizabilities at the magic wavelengths. The resonant wavelengths Arcs for transitions contributing 
to the npj ac polarizabilities and the corresponding absolute values of the electric-dipole matrix elements are also listed. The 
wavelengths (in vacuum) are given in nm and electric-dipole matrix elements and polarizabilities are given in atomic units. 



Rb 


Transition: 5p3/2mj 


- 5s 






= 1/2 




= 3/2 




Resonances 


d 


'^rcs 


-^magic 


'^('^magic ) 


-^magic 


( Aniagic ) 




5P3/2 - 4d5/2 


10.90(9) 


1529.37 












5p3/2 - 4d3/2 


3.63(3) 


1529.26 


1414.8(5) 


456(1) 








5P3/2 - 6s 


6.05(3) 


1366.87 


792.00(1) 


-6910(30) 


789.98(2) 


125(35) 




5p3/2 - 5s 


5.977(4) 


780.24 












5P3/2 - 5d3/2 


0.67(3) 


776.16 












J"5/2 


1.98(9) 


775.98 


775.84(1) 


-19990(70) 


775.77(3) 


-19700(130) 




^Tl.^ m 7 Q 

'-*P3/2 ' 


J- .OiJU \^ J 


741.02 


740.07(1) 


-2494(4) 








^P3/2 '-*'*3/2 


0.51(5) 


630.10 












5P3/2 - 6d5/2 


1.51(15) 


630.01 


697 Sf'i'l 


-533(4) 


696 9(Q'\ 






5P3/2 - 8s 


0.708(2) 


616.13 


614.7(1) 


-477('1 ") 






Cs 


Transition: 6p3/2mj 


- 6s 






= 1/2 


\iTh\ 


= 3/2 




6p3/2 - 5d3/2 


3.166(16) 


3614.09 


3611.8(2) 


422(1) 


3589(1) 


422(1) 




6P3/2 - 5d5/2 


9.59(8) 


3490.97 


1910(6) 


498(2) 








6P3/2 - 7s 


6.47(3) 


1469.89 


932.4(8) 


3197(50) 


940.2(1.7) 


2810(70) 




6P3/2 - 6d3/2 


2.10(5) 


921.11 


921.01(3) 


4088(10) 


920.18(6) 


4180(14) 




6P3/2 - 6d5/2 


6.15(14) 


917.48 


887.95(10) 


-5600(100) 


883.4(2) 


-1550(90) 




6P3/2 - 6s 


6.324(7) 


852.35 












6P3/2 - 8s 


1.461 


794.61 
















793.07(2) 


-2074(5) 








6P3/2 - 7d3/2 


0.976(9) 


698.54 


698.524(2) 


-697(2) 


698.346(4) 


-696(2) 




OP3/2 '"5/2 




fiQ7 '\') 


687.3(3) 


-635(3) 


684.1(5) 


-618(4) 




6p3/2 - 9s 


0.770 


658.83 


657.05(9) 


-500(1) 








6P3/2 - 8d3/2 


0.607(8) 


621.93 


621.924(2) 


-388(1) 


621.844(3) 


-388(1) 




6P3/2 - 8d5/2 


1.81(2) 


621.48 


615.5(8) 


-371(3) 


614(3) 


-367(8) 




6p3/2 - lOs 


0.509 


603.58 


602.6(4) 


-339(1) 







C. Rb 

We list the magic wavelengths Amagic above 600 nm for 
the 5pi/2 — 5s transition in Rb and the 6pi/2 — 6s tran- 
sition in Cs in Table fVl In this case, all Rb 5pi/2 — nlj 
resonances have significant spacing allowing us to deter- 
mine the corresponding magic wavelengths. The magic 
wavelengths above 600 nm for the 6^3/2 — 5s transition in 
Rb and 6^3/ 2 — 6s transition in Cs are grouped together 
in Table IVll The transitions up to 5p — 10s and bp — 9d 



are included in the main term calculation of the Rb 5p 
polarizabilities and the remainder is evaluated in the DF 
approximation. The bp — 5s matrix elements are taken 
from Ref. [3, and the bp— Ad El matrix elements are the 
recommended values derived from the Stark shift mea- 
surements in Ref. As we discussed in Section [Til 
the correlation correction is very large for the bp — bd 
transitions; the DF values for the 5^3/2 — 5^5/2 transi- 
tion is 0.5 a.u. while our final value is 2.0 a.u. However, 
nearly entire correlation correction to this value comes 
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FIG. 6: ac Stark shifts for the transition from 5p3/2F' = 3M' 
sub levels to 5sFM sub levels in Rb as a function of wave- 
length. The electric field intensity is taken to be 1 MW/cm^. 
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FIG. 7: Polarizability of the 6P3/2 state and polarizability 
of ground state of Cs as a function of wavelength. Magic 
wavelengths for the 6P3/2 — 6s transition of Cs are found to 
be at 932 and 940 nm depending on the m-, value. 



from the single all-order term which can be more ac- 
curately estimated by the scaling procedure described in 
Refs. ^, 1I&, 16]. To evaluate the uncertainty of these val- 
ues, we also conducted another calculation including the 
triple excitations relevant to the correction of the domi- 
nant correlation term (SDpT method) , and repeated the 
scaling procedure for the SDpT calculation. We took the 
spread of the final values and the SDpT ab initio and 
SDpT scaled values to be the uncertainty of the final 
numbers. Nevertheless, even such an elaborate calcula- 
tion still gives an estimated uncertainty of 4.5%. 

We illustrate the Amagic for the 5p3/2 ^ 5s transition 
near 791 nm in Fig. [5l We note that this case is different 
from that of Na illustrated in Fig. [31 where both ao -f- 
a2 and ao — a2 curves for the 8^3/2 polarizability cross 
the 3s polarizability curve at approximately the same 
polarizability values. In the Rb case near 791 nm, ao + 
a2 and ao — a2 curves for the 5p3/2 polarizability cross 
the 5s polarizability curve at 125 a.u. and -6910 a.u., 
respectively. As a result, the \M'\ = 3 curve on the ac 
Stark shift plot for the transition between hyperfine sub 
levels shown in Fig. [6] is significantly split from the curves 
for the other sublevels. 

The magic wavelengths for the 5p3/2 — 5s transition be- 
tween the fine-structure components of the 5p—ndj levels 
are not listed owing to very small fine structures of these 
levels. We note that crossings for all rrij sublevels should 
be present between the fine-structure components of the 
bp — ndj lines. We illustrate such magic wavelengths 
for Cs, which has substantially larger ndj fine-structure 
splittings. 



D. Cs 

Our results for Cs are listed in Tables |V] and IVII The 
values of the 6p — 6s matrix elements are taken from (35i] , 
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FIG. 8: ac Stark shifts for the transition from 6P3/2F' = 5M' 
sub levels to 6sFM sub levels in Cs as a function of wave- 
length. The electric field intensity is taken to be 1 MW/cm'^. 

and the values for 6p — 7s transitions are taken from 
the results compiled in [l^ (derived from the 7s life- 
time value). We derived the 6pi/2 — 5d3/2 value from 
the experimental value of the Dl line Stark shift in Cs 
^\ combined with the experimental ground state polar- 
izability value from [soj . The procedure for deriving the 
matrix element values from the Stark shifts is described 
in Ref . [22| . We use the theoretical values of the ratios of 
the 6pi/2 - 5(^3/2, 6p3/2 - 5(^3/2, and 6^3/2 - 5dr,/2 values 
from Ref. [1] to obtain the values for the 6^3/2 — 5^3/2 and 
6p3/2 — 5^5/2 matrix elements. We use the experimental 
energy levels from [10, 11, 36j, and references therein to 
obtain the resonance wavelength values. The transitions 
up to 6p — 12s and 6p — 9d are included into the main 
term calculation of the polarizabilities and the remainder 
is evaluated in the DF approximation. 

We find that there are no magic wavelengths for 
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the Qpi/2 — 6s transition in between the Qpi/2 — 6s, 
6pi/2 — 6(i3/2, and 6pi/2 — 8s resonances whereas there 
are the magic wavelengths in between the corresponding 
resonances in Rb. The difference between the Rb and Cs 
cases is in the placement of the 6^3/2 — 6s resonance in 
Cs and 5p3/2 — 5s resonance in Rb. In Rb, 5p3/2 — 5s 
resonance is at 780 nm and follows the 5pi/2 — 5s one. In 
Cs, the 6p3/2 — 6s resonance is at 852 nm and is located 
in between the &P1/2 — 6^3/2 and 6pi/2 — 8s resonances 
owing to much larger Qp fine-structure splitting. As a 
result, there are no magic wavelengths in this range. 

Also unlike the Rb case, the magic wavelengths around 
935 nm for the 6^3/2 — 6s transition in Cs correspond to 
similar values of the polarizability (about 3000 a.u.) for 
all sublevels as illustrated in Fig. [T] The nearest reso- 
nances to this magic wavelength are 6^3/2 — 6dj ones; 
therefore the contributions from these transitions are 
dominant. To improve the accuracy of these values, we 
conducted a more accurate calculation for these transi- 
tions following the hp — 5d Rb calculation described in 
the previous subsection. As a result, we expect our val- 
ues of the 6p — 6d matrix elements to be more accurate 
than the one quoted in Ref. Nevertheless, the uncer- 
tainties in the values of the corresponding magic wave- 
lengths are quite high because the 6s and 6^3/2 polariz- 
ability curves cross at very small angles. As a result, even 
relatively small uncertainties in the values of the polar- 
izabilities propagate into significant uncertainties in the 



values of the magic wavelengths. Our values for these 
magic wavelengths are in good agreement with previous 
studies [1,[33. The ac Stark shift of the 6^3/2-^' = 5M' 
to 6sFM transition as a function of wavelength at the 
925 — 945 nm range is plotted in Fig. [H 



V. CONCLUSION 

We have calculated the ground ns state and np state ac 
polarizabilities in Na, K, Rb, and Cs using the relativis- 
tic all-order method and evaluated the uncertainties of 
these values. The static polarizability values were found 
to be in excellent agreement with previous experimental 
and theoretical results. We have used our calculations to 
identify the magic wavelengths at which the ac polariz- 
abilities of the alkali-metal atoms in the ground state are 
equal to the ac polarizabilities in the excited npj states 
facilitating state-insensitive cooling and trapping. 
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